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Topological insulators supporting non-abelian anyonic excitations are at the center of attention 
as candidates for topological quantum computation. In this paper, we analyze the ground-state 
properties of disordered non-abelian anyonic chains. The resemblance of fusion rules of non-abelian 
anyons and real space decimation strongly suggests that disordered chains of such anyons generi- 
cally exhibit infinite-randomness phases. Concentrating on the disordered golden chain model with 
nearest-neighbor coupling, we show that Fibonacci anyons with the fusion rule r ® r = 1 © r exhibit 
two infinite-randomness phases: a random-singlet phase when all bonds prefer the trivial fusion 
channel, and a mixed phase which occurs whenever a finite density of bonds prefers the r fusion 
channel. Real space RG analysis shows that the random-singlet fixed point is unstable to the mixed 
fixed point. By analyzing the entanglement entropy of the mixed phase, we find its effective central 
charge, and find that it increases along the RG flow from the random singlet point, thus ruling out 
a c-theorem for the effective central charge. 
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I. INTRODUCTION 



One of the most interesting frontiers of physics is the 
behavior of interacting, many body, quantum systems. 
Such systems are particularly challenging and rich when 
considered in low dimensionalities, and in the presence 
of disorder. A common platform for the discussion of 
collective behavior is a quantum magnet. Already in 
one dimension, where the behavior of quantum magnets 
is supposed to be the simplest, surprises emerged: the 
Haldane gap of integer spin chains^'^, and, more impor- 
tantly for this work, the random singlet phase^'"'. The 
latter describes the ground state of a disordered spin-1/2 
Heisenberg chain, where the spins pair up in and form 
singlets in a random fashion (fig. 1). Most of the singlets 
connect near neighbors, but some are very long ranged, 
and lead to algebraically decaying average correlations. 
The random singlet phase is the first known example of 
the infinite-randomness paradigm of Id random systems. 
Contrary to the quantum scaling in pure systems, where 
1/E ~ L^, infinite randomness systems obey the scaling 
ln_E| ^ L'^ , and exhibit many other intriguing proper- 
ties. 

Another paradigm for interacting quantum matter is 
the fractional quantum Hall system. In addition to ro- 
bust fractionally-charged excitations. Hall bars with elec- 
tronic densities tuned to special fractions, such as i' = 
5/2, and v = 12/5 are expected to exhibit non-abelian 
quasi-particles excitations and defects^ '^, which may be 
used to reahze a topologically protected qubit^'*, but 
more importantly, provide an example of a completely 



new type of quantum matter. Non-abelian anyons dis- 
play a remarkable feature: the dimension of the Hilbert 
space spanned by N non-abelions grows asymptotically 
as D^ , where D, the quantum dimension of the non- 
abelion, is irrational. This is a consequence of the so- 
called fusion rules of the non-abelions. In this paper we 
will investigate the properties of one-dimensional disor- 
dered systems composed of non-abelions, concentrating 
on the case of the Fibonacci anyons, for which the al- 
lowed values of the so-called topological charge can be 
either 1 or t. 

The investigation of the so-called Fibonacci chain has 
so far concentrated on the analysis and phases in the 
translationally invariant case^'^". It turns out that the 
system is exactly solvable by mapping to an RSOS model, 
and is described at low energies by a minimal model con- 
formal field theory with central charge c = ^^ in the 
antiferromagnetic case (favoring fusion into the trivial 
channel) and c = ^ in the ferromagnetic case (favoring 
fusion into the r channel) . The richness of this example 
stems from the unique structure of the Hilbert space of a 
system comprising non-abelian anyons. As we argue be- 
low, however, an important insight is that the construc- 
tion of the Hilbert state of a random non-abelian chain 
is analogous to the construction of the ground state and 
low-lying excitation spectrum of a spin-chain. Further- 
more, contrary to spin chains comprising garden-variety 
spins, the Hilbert space structure of non-abelian chains 
guarantees the appearance of an infinite randomness scal- 
ing in the presence of disorder. The exotic nature of 
the non-abelions suggests that the infinite-randomness 
phases they will exhibit will be new, and perhaps even 
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FIG. 1: Sample random singlet ground state of a spin-1/2 an- 
tiferromagentic chain. Each line represents a singlet. Singlets 
form in a random fashion, mostly between nearest neighbors, 
but occasionally between distant sites giving rise to an aver- 
age correlation that decays algebraically. 



expand our dictionary of infinite-randomness universality 

classes, currently limited to the permutation symmetric 

1 1-1 s 
sequence ". 

Another interesting aspect of infinite-randomness 
phases is their entanglement entropy. The bipartite en- 
tanglement entropy of a pure spin chain at criticality 
scales logarithmically with its size, and is proportional 
to the central charge of the conformal field theory de- 
scribing the critical point ^^"^^. Furthermore, the central 
charge, and therefore also the entanglement entropy of 
a pure spin chain, obeys the Zamolodchikov c-theorem: 
it must decrease along renormalization group fiow lines. 
Random spin chains also have an entropy that scales log- 
arithmically with size, and with a universal coefficient 
that we identify as an effective central charge^^"^°. An 
outstanding question has been whether the c-theorem ap- 
plies to renormalization group flows between infinite ran- 
domness fixed points of the random chains. The evidence 
so far has been limited, since the overwhelming major- 
ity of entanglement entropy calculations were done in the 
random singlet phase of various systems. The only ex- 
ception so far has been the entanglement entropy at the 
critical point between the Haldane phase and the ran- 
dom singlet phase of a spin-1 random antiferromagnet, 
where the effective central charge indeed decreases along 
real-space RG fiow lines^'*. In this paper, we find a real- 
space RG flow along which the effective central charge 
increases, thus violating any conjectured c-theorem for 
flows between strong randomness fixed points. 

The first analysis of the random Fibonacci chain^^ 
concentrated on the random antiferromagnetic Fibonacci 
chain. A random singlet phase was found with an effec- 
tive central charge reflecting the quantum dimension of 
the Fibonacci anyons: Cgff = Inr with r — ^^^ being 
the golden ratio. In this paper we extend the analysis to 
the completely disordered Fibonacci chain, which con- 
tains both AFM (antiferromagnetic, favoring fusion into 
a singlet) as well as FM (ferromagnetic, favoring fusion 
into a r) bonds - we will refer to it as the mixed Fibonacci 
chain. We flnd that the random singlet phase is unstable 
to FM perturbations, and flows to a stable flxed point 
which, at low energy, is described by an equal mixture of 
FM and AFM bonds, with identical inflnite randomness 
universal strength distributions. We calculate the entan- 
glement entropy of this new flxed point, and show that it 
is larger than that in the unstable AFM random-singlet 
point. We thus have an example of a flow between two 
inflnite randomness flxed points along which the effective 



central charge increases. 

In the rest of the paper we will describe our analysis 
of the random Fibonacci chain. In Sec. II we review the 
real-space renormalization group and the Hilbert space 
and Hamiltonian of the Fibonacci chain. In Sec. Ill we 
define the model, and discuss the decimation rules neces- 
sary for a real-space RG analysis. The decimation steps 
are then used for the calculation of flow equations for the 
disorder distribution. We use these to derive the phase 
diagram and investigate the stability of the flxed points 
found from the flow equations. Sec. IV will describe 
the entanglement entropy calculation for the random Fi- 
bonacci chain. Before concluding, we will discuss the 
correspondence between the construction of the ground 
state of a random spin chain, and the construction of the 
Hilbert space of a chain of non-abelian anyons. This pro- 
vides the basis for further investigation of other kinds of 
non-abelian chains, such as the full SU{2)k sequence. 



II. BACKGROUND 
A. Real Space Renormalization Group 

To flnd the ground state of disordered spin chains, 
Ma and Dasgupta introduced the strong disorder real- 
space renormalization group^^'^'^. The random spin 1/2 
Heisenberg model provides the simplest example for this 
method. The model is given by: 



H — 2_^ Jis+i Si • S.i-|_i, 



(1) 



where the couplings Ji^i+i > are positive and randomly 
distributed. Note that, as far as the Hilbert space is 
concerned, we have for two neighboring sites 



1 



= 0©1, 



(2) 



and the interactions in the Hamiltonian simply give an 
energy splitting between the two representations on the 
right hand side. The procedure now is to pick the largest 
Ji,i+i, which effectively truncates the excited triplet and 
leaves the ground state in a singlet, and do perturbation 
theory around that state. Quantum fluctuations then 
induce an effective coupling according to the so-called 
Ma-Dasgupta rule^^^^^: 



J^-l 



Ji-l,iJi+l 



+ 1,1+2 



2J, 



i,i+l 



(3) 



So sites i and i + 1 are decimated and replaced with an 
effective interaction between i — \ and i + 2. Iteration of 
this procedure produces bonds on all length scales. This 
is the random singlet ground state. 

A quantitative description is obtained by tracking the 
RG flow of the coupling distribution. It is useful to em- 
ploy logarithmic couplings'^: 

A,.+i = hi (4) 



J, 



i,i+l 



where il = max^Ji^i+i. 
Dasgupta rule (3) reads 



In these variables the Ma- 



a) 



Pi-l,i+2 — /3i-l,i + /3i+ 



1,1+2 



(5) 



(up to an additive constant of In 2 which can be safely 
neglected). As the couplings get decimated V, decreases. 
It is convenient to define the RG flow parameter as 



r = ln 



n 



(6) 



where VIq is the maximal coupling of the bare Hamilto- 
nian. Let Pr{P) be the distribution of couplings. We can 
derive a flow equation for Pr(/3) by decimating the cou- 
plings in the infinitesimal interval f3 = [0,dT] and seeing 
how their probabilistic weight is redistributed. We obtain 



"^ P(B)- ^^^ 



(7) 



/•OO /'OO 

P(0) / d(3j d/32<5(/3-/3i-/?)Pr(/3i)Pr(/32) 
Jo Jo 

The first term comes from the overall change of scale, and 
the second from the Ma-Dasgupta rule. These equations 
have a solution 



PrW = ^e-'^/^ 



(8) 



which is an attractive fixed point to essentially all phys- 
ical initial configurations. This solution permits us to 
read off features of the random singlet phase; for example 
one can with a little more work derive the energy-length 
scaling relation: 



Li/2_r = ln(l/P). 
which thus has the exponent: 

V- = 1/2. 



(9) 



(10) 



B. Hilbert Space and Hamiltonian of the Fibonacci 

Chain 



We now construct the Hilbert space and Hamiltonian 
of the Fibonacci chain. The system is modeled as a chain 
of non-abelian anyons carrying the non-trivial topological 
charge r. Heuristically, we want the property that 



> T — 1 ' 



(11) 



which states that the Hilbert space of two neighboring 
r's is the direct sum of a trivial component and another 
copy of T. This unusual property immediately prevents 
us from describing the Hilbert space as a tensor product 
of local degrees of freedom. Indeed, a naive interpreta- 
tion of the tensor product in (11) would give the dimen- 
sion of the space r to be the golden ratio, an irrational 
number. This problem is resolved by the adoption of the 
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FIG. 2: a) A state in the Hilbert space: the labels 1 and 
r specify the total topological charge of all the sites to the 
left (or equivalently right) of the bond, with the fusion rules 
obeyed at each trivalent node b) The same state in graph 
notation - we only draw the r's. 



machinery of truncated tensor products of representations 
of SU{2) at level k, but rather than developing it here we 
instead give two elementary constructions of the Hilbert 
space. We note, however, the analogy between (2) and 
(11); indeed the Hamiltonian, defined below, will simply 
yield an energy splitting between the two representations 
on the right hand side of (11). 

The simplest way to construct the Hilbert space is to 
define basis states by labeling each link between two r's 
with a 1 or r, with the constraint that one is not allowed 
to have two consecutive I's- fig. 2(a). The dimension Z? at 
of the Hilbert space for N sites then follows the Fibonacci 
recursion 



Dn = D 



w-i 



D 



N-2 



(12) 



which is solved by Dn ^ r . Thus there are r ~ 1.618 
"degrees of freedom" on each site (note: we use r to de- 
note the nontrivial topological charge and the value of the 
golden mean, as well as the corresponding representation 
of SU{2)k where appropriate). 

While this link description of the Hilbert space is most 
convenient computationally, there is an equivalent but 
more abstract one that is useful in defining the Hamilto- 
nian and carrying out the real space RG procedure. In 
this abstract description the Hilbert space is defined as 
the set of all trivalent graphs with endpoints at the N 
nodes, modulo the P-matrix relations - see fig. 3 - where 
the F matrix is 



P: 



-1/2- 



(13) 



The edges of the graph represent nontrivial topological 
charge r and trivalent vertices represent the fusion of two 
r's into another r. To relate this graphical picture to the 
link basis, note that the link basis states can be viewed as 
trivalent graphs, as in fig. 2(b), and any other trivalent 
graph can be reduced to a superposition of these using 
P-matrix moves (for example, see fig. 4). The inner 
product of two graphs is defined by refiecting one of the 
graphs and concatenating it with the other along the N 
nodes. 
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FIG. 3: The F-matrix relations, a) graphs that can be dis- 
connected by cutting one edge are equal to (the no tadpole 
condition) and disconnected loops are worth r. b) local re- 
connection rules are given by the .F-inatrix, defined in the 
text. Here m and n are binary variables equal to either r or 
1 - i.e., the link is either there or not 
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FIG. 4: The trivalent graph on the left hand side expressed 
in terms of the link basis. 
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FIG. 5: a) Projection operator P/* on a pair of sites, b) 
Projection operator Pf on a different pair of sites. 




FIG. 6: A rainbow diagram, 
middle. 



The partition bond is in the 



We now define the Hamiltonian. There are two kinds 
of nearest neighbor interactions we consider: we can ei- 
ther project onto total topological charge r of the pair, 
in which case we refer to the interaction as ferromagnetic 
(F), or onto the trivial charge 1 - this interaction is an- 
tiferromagnetic (A). The F and A designations are by 
analogy with the spin 1/2 case, where antiferromagnetic 
interactions favor a singlet, which has trivial spin, and 
ferromagnetic interactions favor non-zero total spin 



n = J2Mi-pl 



(14) 



where for each site Ji is a positive random number with 
a given distribution, and Jli — A for the Hamiltonian 
describing the AFM fixed point, while Y^i = F, A at ran- 
dom for the Hamiltonian describing the mixed FM/AFM 

fixed point. Here Pj are the projectors onto the fer- 
romagnetic and antiferromagnetic sectors of the pairs of 
sites. These projection operators are defined graphically 
in fig. 5. They can also be viewed in the link basis if 
we apply _F-matrix rules to decompose the result of the 
concatenations in fig. 5. 

As an illustration of the complicated nature of the 
Hilbert space, we analyze the "rainbow" state that arises 
in the entanglement entropy calculation for the AFM 
chain^^ - see fig. 6 (a). If we assume the partition bond 
(i.e., the bond that divides the system into two subsys- 
tems) to lie in the middle of the chain, then this state 



turns out to be quite entangled. We will be more precise 
later, but the idea is that to compute the entanglement 
entropy, we want to use P-matrix moves to write the 
state as a superposition of states shown in fig. 6 (b). 
This way we push all the nontrivial parts of the graph 
into one of the two halves, and the entanglement can be 
read off from the coefficients of the new states. Carrying 
this out, Ref. 21 showed that for large number TV of sin- 
glets, the entropy is asymptotically N log2 r. Thus the 
asymptotic contribution of each singlet is log2 r. Note, 
however, that when N is small, there are deviations from 
this form. In particular, when A^ = 1, so that we have 
only one t anyon in each half, the Hilbert space has only 
one state and so the state counting entropy is 0. These 
sorts of subtleties will be treated carefully when we do 
our entropy calculation for the mixed FM/AFM fixed 
point. 



III. INFINITE RANDOMNESS FIXED POINTS 
OF THE GOLDEN CHAINS 

A. Ma-Dasgupta decimation rules 

The graphical description of the Hilbert space in terms 
of trivalent graphs shows how the real space RG method 
can be generalized to the case of the mixed Fibonacci 
chain, containing both FM and AFM interactions. As 



before one first picks the largest coupling Ji in eqn. (14) 
and assumes that it localizes a state on i and i + l with to- 
tal topological charge either 1 or r, depending on whether 
the interaction is AFM or FM. Graphically this localiza- 
tion is just a restriction to graphs that have a singlet 
spanning the two sites (AFM case) or graphs that have 
the two r's at i and z -I- 1 fuse into another r (FM case). 
Again the state of the two sites is perturbed by the other 
two bonds connecting these sites to the rest of the chain. 
To study the effect of this perturbation consider the 
Hamiltonian (14) acting on four tau particles with site 
labels 1 through 4. Using the fact that Pf = 1 — Pf^ this 
Hamiltonian can, up to an irrelevant constant, be taken 
to be 



H = -JiP^ - J^P^ - .hP^, 



(15) 



where now the sign of a given Ji, connecting particles at 
sites i and i+l, can be positive or negative, corresponding 
to AFM or FM bonds, respectively. We then assume that 
J2 is the highest energy bond, with | J2 1 ^ | >/i | , | '/2 1 , and 
write (15) as H = Hq + H' where Hq = -J2P2 is the 



-JiP^ 



,hP^ is 



^3^3 



"unperturbed" Hamiltonian and H' 
the perturbation. 

First consider the case of decimating an AFM bond for 
which J2 > 0. The two degenerate ground states of Hq 
will have a singlet connecting particles 2 and 3 (i.e., these 
two particles will have total topological charge 1) while 
particles 1 and 4 can combine to either have topological 
charge 1 or r. We denote these two unperturbed states 
lipi) and IV't)- Since the total topological charge of the 
four particles is a "good" quantum number, the pertur- 
bation H' will lead to an energy splitting, Je// , between 
the state of these four particles with total topological 
charge 1 and the state with total topological charge t. 
This energy splitting can then be described by a new 
Hamiltonian H^ff — —JeffPi, where now P^ acts on 
particles 1 and 4. It is straightforward to compute Jeff 
using second order perturbation theory with the result, 

Jeff = 7 7 ■ (16) 



J2 



J2 



In this expression the FM projection operator P^ 
projects TJ'IV'i.t) onto the excited Hilbert space of the 
unperturbed Hamiltonian Hq with energy J2 above the 
ground state. 

Using the techniques described in the previous section 
the matrix elements appearing in (16) can be evaluated 
to find 



lii^^lH'Pi^H'li^,)]' = {Ji + Jsf^ 



- (J1 + J3) 



1 



\{i^r\H'PiH'\ijr)\' = {Jt + JiH 



(Ji + Js) — 



(17) 



(18) 
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FIG. 7: Particles 2 and 3 are decimated ferromagnetically. 
The efltect on an AFM projection acting on 1 and 2 is simply 
to turn it into an FM projection. For an AFM projector (FM 
projector) there is a numerical factor associated to the graphi- 
cal representation above of r^^ (r^^"). Combining these, we 
get r~^(r~^")^ — t~^ , of which r~^" gets absorbed in the 
normalization for the new FM projector and another r~^" 
goes into the normalization of the wavefunction, leaving r~^. 
Thus P on sites 1 and 2 turns into -P^ . A similar graphical 
argument shows that P^ turns into -P' . 



It then follows that 



J< 



eff 



2 J1J3 

t2 J, 



(19) 



Thus we see that when an AFM bond is decimated the 
usual Ma-Dasgupta rule holds. The value of the coef- 
ficient 2/t^ is not significant except for the fact that, 
because it is less than 1, Je// will always be less than J2. 
As for the usual Ma-Dasgupta rule, the resulting inter- 
action will be AFM if Ji and J3 have the same sign, and 
FM if Ji and J3 have opposite signs. 

Next consider the case of a FM bond for which J2 < 
in (15). In this case the two tau particles connected by J2 
fuse to form a cluster with topological charge r. The tau 
particles on either side of this cluster will then interact 
with it, but with modified interaction strengths Ji and 

J3. 

To compute these modified interactions, consider par- 
ticle 1, which is coupled to the newly formed cluster 
through the "bare" interaction Ji. To see the effect of 
the decimation on any operator O on the pair of parti- 
cles 1 and 2, we simply project this operator down to the 
decimated subspace: Onow = P^ O P^ where P^ acts on 
particles 2 and 3. Composing the operators graphically 
in fig. 7 we see that P^ turns into P^ and vice versa, 
with an extra factor of 1/r. 

Another way to see this is as follows. There are two 
possible states for the tau particle at site 1 and the r 
cluster formed by particles 2 and 3 — the total topo- 
logical charge of all three particles can be either 1 or 
T. Again we denote these two states |V'i) and IV't)- In 
this case the effective interaction can be computed using 
first order perturbation theory in H' — ^JiP^ with the 
result 

Jl = {lPi\JiP{'\^Pi) - (V'rlJlPi^lVr). (20) 

The calculation of these matrix elements is again 



straightforward and wc find that 

(V'llPi^lVi) = 0, 



(V'rlPi^lV'r) 



1 



Thus we obtain 



•h 



-Ji. 



(21) 
(22) 

(23) 



The same argument imphes that J3 — — - J3. The essen- 
tial feature here is that when a FM bond is decimated 
the sign of the effective interaction with the neighboring 
tau particles is flipped — FM bonds become AFM bonds 
and vice versa. In addition there is numerical reduction 
of the bond strength by a factor of 1/t. However, as 
above, the value of this coefficient is not important for 
determining the fixed-point behavior of the model — the 
only important fact is that it is less than 1 so we are 
guaranteed that | Ji| < | J2I. 

The two above results for decimation of strong bonds 
constitute the strong-randomness RG rules of the mixed 
random Fibonacci chain. In the following we will derive 
the RG flow equations for this case and show that it has 
non-trivial fixed points. 



B. Flow equations for the Fibonacci chain 

In order to explore the phase diagram of the golden 
chain, we must first turn the Ma-Dasgupta rules for the 
decimation of ferromagnetic (FM) and antifcrromagnetic 
(AFM) bonds, eq. (19), to fiow equations. This goal was 
partially achieved in Ref. 21 for a golden chain which 
contains only AFM bonds. As we shall see, including FM 
bonds in this analysis reveals a new fixed point, where 
the number of FM and AFM bonds is the same. 

We begin our analysis by introducing the logarithmic 



J 



notation for bond strengths: 



A = In ■ 






(24) 



where J7 = maxi{|Ji|} so that the AFM Ma-Dasgupta 
rule (19) reads (3eff = Pi-i +/3i+i — In C. Note that while 
the /?i's carry the information about bond strengths, they 
do not specify whether a bond is FM and AFM, and 
(3i > 0. Let us next define the coupling distributions for 
the AFM (positive J's) and FM bonds (negative J's), 
respectively: 



P(/?), N{(3). 



(25) 



The probability of a bond to be AFM (A), or FM (F) are 
thus 







(26) 



In addition, we define F = In -jy- to be the logarithmic 
RG fiow parameter. Its initial value is a non-universal 
constant of order 1, and as the RG progresses, it flows to 
00. 

The flow equations for P{(3) and N{l3) are derived in 
analogous fashion to those of the distributions in the spin 
1/2 problem'^. Roughly speaking, the terms appearing in 
the two flow equations are the result of: (a) rescaling 
of the UV cutoff ri, (b) decimation of an AFM bond, 
(c) decimation of a FM bond. Below we will write the 
fiow equations with each term followed by an explana- 
tion or a diagram of the decimation step giving rise to 
it. In the graphical representation on the right column 
below, A and F represent antifcrromagnetic and ferro- 
magnetic bonds respectively, and the bond decimated is 
represented by the bold letter with the hat. Let us start 
with the flow of the AFM bond distribution: 



dP 

dV 



dP 

dfi 



CXD 



00 



+P(0) / dPi J dp25{fil + /32 - (i)P{Pl)P{p2) 





00 00 



-P(0) / d/3i / d/32(5(/3i + (i2- l3)N{h)N{ti2) 


-2P(0)P 

+2N{0)N 

-2N{0)P 

+ {2P{0)+N{0))P 



(cutoff rescaling) 

•P»A»P»^»^» 

(neighbor removal in AFM decimation) 

• F • F • F»^» A • A» 

(removal of neighboring AFM in FM decimation) 



(27) 



The last term feeds back the probability of bonds lost due to an AFM decimation, which removes a net of two bonds 
[2P(0)], and due to a FM decimation, which removes a single bond [A^(0)]. 



Carrying out the analogous considerations for the FM bond distribution: 

^ = ^ (cutoff rescaling) 

+2P(0) J dpi Jd(32S{(3i +P2- P)N{l3i)P{p2) .i^.A.A.^.F. 



-2P{0)N 

+2N{0)P 

-~2N{0)N 

-{2P{0) + N{0))N. 



(neighbor removal in AFM decimation) 

• A»F»A»^»F • F • 

(removal of neighboring FM in FM decimation) 



(28) 



Once more, the last line makes sure that probability is conserved. 

Adding up all the above terms yields the following concise flow equations: 



dP _ dP 



^ + P(0) {P (g> P + N (g, N)) + 2N{0)N{I3) - 7V(0)P(/3) 



^ = ^ + 2P{0)N (g>P- N{0)N{I3) + 2N{0)P{I3), 
where we also introduce the notation: 



(29) 



F®G^ I dxi / dx25{x - xi - X2)F{xi)G{x2). 


I 



(30) 



C. Fixed points of the real-space RG 

From the flow equations, Eqs. (29), we can find the 
fixed points of the golden chain. These appear as attrac- 
tors of the integro-differential equations. To find them, 
we first note that we can eliminate F by guessing a scale- 
invariant solution"^: 

Pt{P) = \p{P/n NriP) - ^n(/3/F). (31) 

Substituting this scaling ansatz gives: 

—p={l + x)p' + po(p <E) p + n(E) n) + 2non — nop 
—n = (1 + x)n' + 2pQn ® p -\- 2nop — non. 

(32) 
Furthermore, the convolution hidden by the ® sign com- 
pells us to assume an exponential form for the unknown 
functions n(x), p{x): 



p{x) ~ Pqc '^^ n(x) = hqc '^^ . 



This ansatz reduces the integro-differential equations, 
Eqs. (29), to a set of three simple algebraic equations: 



7=pl + nl 
no(7-2p2)=0 



(34) 



The exponential ansatz and the resulting Eqs. (34) 
reveal two fixed-point solutions. A first solution of (34) 
corresponds to the pure AFM fixed point: 



no 



0. po = 7 = 1- 



(35) 



This is the random singlet phase discussed in Ref. 21. A 
new fixed point, however, is found by allowing no to be 
nonzero: 



7 = 2, po = no = 1. 



(36) 



This fixed point has an equal proportion of FM and AFM 
bonds, and although it is an infinite-randomness fixed 
point, it is not a random-singlet point. Translating back 
to the original variables, the coupling distributions are: 

7V(/3)=P(/3) = ie-2'3/r^ (37) 



While 7 = 2 is one universal critical exponent de- 
scribing the universality class of the mixed FM/AFM 
phase, another critical exponent is -0, which describes 
the energy-length scaling: 



ml 

E 



L'f'. 



(33) This is equivalent to: 



pi/^- 



(38) 



(39) 



with n ^ 1/L here being the density of undecimated 
sites. 

To obtain ip, let us compute the density of free sites 
at the energy scale F. A FM bond decimation eliminates 
one site, while an AFM decimation eliminates two sites. 
This implies that the total density of undecimated sites 
obeys: 



1 dn 
n dV 



-N{0) - 2P(0) 



1 



(40) 



and therefore: 



n = 



ni_ 



(41) 



which corresponds to the infinite randomness critical ex- 
ponent: 



^ = 1/3. 



(42) 



From ip and 7 of the mixed FM/AFM fixed point of 
the Fibonacci random chain, we see that it is in the same 
universahty class as the fixed point separating the gapped 
Haldane-phase and the random singlet phase of the S = 1 
random chain. 



D. Stability of the phases 

In the previous section we found the two fixed points 
of the random golden chain. In order to construct its 
phase diagram, however, we must also study the stabil- 
ity of these fixed points. As it turns out, the random- 
singlet phase is actually unstable^ and flows to the mixed 
FM/AFM fixed point. 

Let us begin our analysis with the mixed FM/AFM 
phase. Assume a perturbation that breaks the balance 
between FM and AFM bonds: 



N{I3) = (1 - 5)ie-2/3/r^ p(/3) ^ (1 + 5)^e-2/5/r. 
Substituting in Eqs. (29) very readily yields: 



(43) 






(44) 



indicating stability with respect to FM/AFM imbalance. 
By establishing the stability of the mixed phase, we es- 
sentially doom the random singlet phase to be unstable. 
Complementing the analysis, however, near the AFM 
random singlet fixed point, we assume: 

N{p)=5^e-f^'^ P{l3) = {l-5)^e-f^/^. (45) 

Again, substitution in Eqs. (29) yield: 



r^^ 9A 



(46) 



which means that FM bonds are a relevant perturbation. 
In addition, we find that the cross-over exponent is: 



X = 2. 



(47) 



These results allow us to draw the fiow diagram. Fig. 
8. Thus the golden chain is in the AFM random-singlet 
phase when it initially consists of only AFM bonds. On 
the other hand, any finite density of FM bonds leads to 
the mixed FM/AFM fixed point. 







-e- 



1/2 



FM/AFM 



FIG. 8: Flow diagram for the mixed Fibonacci chain. Two 
fixed points exist. The AFM fixed point is unstable with 
respect to introduction of FM bonds. A stable fixed point 
exists at the symmetric FM-AFM point. 



IV. ENTANGLEMENT ENTROPY AT THE 
SYMMETRIC FM-AFM POINT 

A. OvervieAV of the calculation 

In this section we calculate the asymptotic scaling of 
the block entanglement entropy of the disordered Fi- 
bonacci chain, that is, the entanglement entropy between 
a region of L consecutive sites and its complement. Be- 
cause of the non-local nature of the Hilbert space, some 
subtleties arise. Let us first define entanglement entropy, 
and then motivate our definition. Given two regions A 
and -B, we have, as illustrated in fig. 9, superselection 
sectors for the topological charge, with the total Hilbert 
space 



H = Ha"<E) Hb° ® Ha^ ® Hb\ 



(48) 



Here the superselection sectors Ha^ and Hb^ can for- 
mally be thought of as n-point disk spaces. Given a state 
tp G H we decompose it according to (48) asip — ip'^ + ^p^. 
Each of these has a Schmidt decomposition 



V'''-E^/^/®^/ 



(49) 



where the states 77/ (g) x/ have unit norm in H. We now 
define the entanglement entropy in the usual way, as 



S = -Y,Xilog,X, 



(50) 



To motivate this definition we note that it is equivalent 
to the standard definition of entanglement entropy when 
we implement the fusion rule constraints via large energy 
penalties E in the Hamiltonian. Specifically, working in 
the link basis for convenience, we enlarge the Hilbert 
space to a space H' that allows all link configurations, 
with terms added to the Hamiltonian to penalize viola- 
tions of the fusion rules. We extend the inner product 
and Hamiltonian to H' in the simplest way possible - say, 
by extending the Hamiltonian to be kE times the iden- 
tity on the space Vk of configurations with k violations of 
the fusion rules, and taking the inner product such that 
Vk is orthogonal to Vi for k ^ I. The new hamiltonian is 
Hermitian on ff', preserves H C H' , and reduces to the 
original Hamiltonian on H . For E much larger than the 
ground state energy in the original problem, the ground 
state and all low lying states in the new problem are the 
same as those in the original one. The upshot is that 
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FIG. 9: Decomposition of the total Hilbert space into super- 
selection sectors. 



H' now has a tensor product decomposition, and entan- 
glement entropy can be defined in the conventional way 
(some care must be taken in normalizing inner products 
on the sub-system Hilbert spaces). This conventional def- 
inition for states in H C H' coincides with ours above. 
We also note that this is how entanglement entropy was 
defined in the numerical algorithm of Ref. 9, which re- 
covered the central charges c = 4/5, 7/10 in the uniform 
case of the golden chain. 

The calculation of the entanglement entropy for the 
Fibonacci chain proceeds along the same lines as pre- 
vious calculations of the entropy. To gain orientation 
for the calculation we will shortly pursue, let us review 
the entanglement entropy calculation for the simplest in- 
stance of an infinite randomness fixed point, the spin 1/2 
Heisenberg model^^. There, the basic idea is to count the 
number of singlets formed over a boundary of the inter- 
val, up to a cutoff size L (energy-length scaling turns this 
into a cutoff in F) . Each singlet contributes 1 to the en- 
tanglement entropy. Real space RG analysis shows that 
the number of singlets is proportional to logF '^ logi, 
so we obtain logarithmic scaling of the block entangle- 
ment entropy. It turns out that this logarithmic scaling 
persists in the disordered fibonacci chain, but obtaining 
the coefficient in front of the log is considerably more dif- 
ficult. For one thing, we already saw earlier (discussion 
preceding fig. 6) that even in the AFM fixed point, ob- 
taining the entanglement entropy required using i<"-moves 
to change to a more convenient basis. 

Obtaining the entropy in the mixed fixed point of the 
Fibonacci chain is even more difficult, because the ground 
state now contains not just singlets but also complicated 
tree-like structures, since two r's can fuse into another r, 
and not just to a singlet. The problem, however, is still 
tractable, although instead of looking at the RG time be- 
tween successive singlets, we must now look at RG times 
between consecutive AFM decimations, and the tree-like 
structures that form between them (fig. 11a). Just as 
each singlet in the AFM case contributed asymptotically 
log2 T in the case of many singlets, we will find a simi- 
lar simplification in the mixed case for a large number of 
AFM decimations - each tree-like structure will asymp- 
totically contribute some amount to the entropy. The RG 
process will average over all trees, so we will have some 
average contribution < S'treo > to the entanglement en- 
tropy. To get the dependence on the block length L, we 
first use energy-length scaling to relate L to the RG flow 
variables: ^\nL ^ — InF ^ /. We show in the next sub- 
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FIG. 10: RG history near the partition bond (represented 
by the dashed vertical line). The partition bond is at some 
point decimated, with equal probabilities of having an AFM 
and FM decimation. In the case of an AFM decimation, we 
get another partition bond, which will at some point be deci- 
mated by an AFM or FM decimation, again with equal proba- 
bility. In the case of an FM decimation, however, we produce 
a partition site, which can decimate, either via AFM or FM 
decimation, with either the site to its left or the site to its 
right. All four of these possibilities are equally likely. Carried 
out further, this process will generate complicated tree-like 
structures. Note that we have ignored decimations that do 
not involve the partition bond or partition site 



section that AFM decimations occur with period 3/2 in 
I. Thus for a block of size L, we have n ^^ ^l ^ ^\nL 
AFM decimations separating tree-like structures strad- 
dling each endpoint of the block, which gives a contribu- 
tion of 



^ <S 
9 < *t- 



> InL. 



(51) 



We will see that there will also be another contribu- 
tion to the entanglement entropy, coming from the resid- 
ual singlets left straddling the endpoints after the tree- 
like structures have been resolved. We will compute this 
" rainbow" contribution carefully later in this section, but 
first we turn to calculating the RG times between the 
various decimations. 



B. RG times betw^een decimations 

A pictorial representation of the RG process is given in 
fig. 10. We see that eventually a ground state of the form 
shown in fig. 11 is generated. To quantitatively under- 
stand real space RG, we will need to compute the (log) 
of the RG times between various types of decimations. 

Now, in the simple case of the spin 1/2 Heisenberg 
model, one is interested in (log of) the RG time between 
successive decimations of the partition bond - the bond 
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through which the boundary of the region passes. Af- 
ter each decimation of the partition bond, the couphng 
distribution at the bond is different than the average dis- 
tribution in the chain. Nevertheless, it is independent of 
the other surviving bonds^'', and in that sense, it 'resets'. 
Thus the RG time duration between successive decima- 
tions obeys a Poisson distribution characterized by one 
number, the average of the log of the RG time between 
successive decimations. No history dependence of this 
number appears for the antiferromagnet. 

Our mixed Fibonacci case is more complicated due to 
the presence of both FM and AFM decimations. Here 
we must consider all possible histories of FM and AFM 
decimations, as illustrated in fig. 10, and compute prob- 
abilities for each. Nevertheless, the Fibonacci chain has 
several simplifying factors which make the problem more 
tractable. First, as in the case of the spin-1 Heiseneberg 
chain calculation, once the partition bond undergoes an 
AFM decimation, the resulting distribution of the cou- 
pling across the partition bond is 'reset', and is indepen- 
dent of the chain's history. Second, the distribution of 
bond strengths is always symmetric with respect to ex- 
change of FM and AFM couplings. It turns out that it is 
characterized by just two numbers, the time between an 
AFM decimation and the next decimation (equally likely 
to be FM and AFM by symmetry), and the time between 
a FM decimation and the next one (again equally likely 
tobeFM and AFM). 

To see this, we notice that the joint probability distri- 
bution of all the bonds takes one of two forms, depend- 
ing on whether we've just had a FM or AFM decimation. 
Immediately following an AFM decimation, we have an 
independent distribution for all bonds, with the partition 
bond having distribution 



Q(/3) 



Y2 



-2/3/r 



(52) 



and all the other bonds following P(/3) = i e-^/^/r. FM 
decimations are even simpler: after an FM decimation, 
all the bonds follow P(/3) ~ jr e^^^^^ . The surrounding 
bonds do get changed from AFM to FM and vice versa, 
but because they are equally likely to be one or the other 
at the mixed fixed point (i.e., the distribution is sym- 
metric with respect to the interchange) there is no net 
effect. One can verify these observations by noting, as in 
17, that under RG evolution following an AFM decima- 
tion the distribution retains its form, with Q changing in 
a complicated way and P evolving as its explicit depen- 
dence on F dictates. In fact, in much the same way as 
is done in Rcf. 17 we can derive an equation for the RG 
evolution of Q: 



We solve it by making the ansatz 



Qr= [a + b 



2f3\ 1 



2/3/r 



r / F 



(54) 



with a and b functions of F. Let I = In^^, where Fn is 

1 

the RG time when the AFM decimation occurred. The 
initial conditions at I — are then a = 0, 6 = 1. Plugging 
the ansatz into (53) then yields 



da 



db 
~dl 



-3a + 26, 



a -26. 



(55) 



The solution is 



-,-1 „-4;\ 



a = -(e--e- 



(56) 



To extract the expected value < I > until the next deci- 
mation, note that the probability p that another decima- 
tion has not occurred by RG time F > Fq is simply 



PT == / d(3Qr{P) =ar + 6r. 



(57) 



The expected value of I at the next decimation is then 



<l> 



I dpi 



2 / aldl 
Jo 



- r{e-'-e-^')dl = -. 



(58) 



Notice that < I > is independent of Fq. 

In a similar manner, we can consider the case where an 
FM decimation of the partition bond has just occurred at 
Fq. In this case all the bond strengths just evolve accord- 
ing to independent distributions P{P) = jr e^^^^^ , so the 
whole situation is characterized by an overall probabilis- 
tic weight w. The RG equation for w is readily derived 
to be 



dw 



-4w 



dV 

F' 



(59) 



where the prefactor of 4 is due to the fact that the site 
containing the partition can be decimated by processes 
on either side of it which can each be either FM or AFM, 
thus leading to four possibilities. This equation is solved 
hy w = (F/Fo)"^ = e^^". The expected value of I at the 
next decimation is then 



<l >= - {de 



.,-4^ 



1 = 4 



^-4:1 



ldl=-. 



(60) 



^^OV(^)-^Qr(^) + 



(53) 



- / dpi dp2 SiP 



Pi-P2)Pr{Pi)Qrm. 



Finally, we compute the average < I > between AFM 
decimations, as follows. First, note that there can be 
any number of FM decimations in between the two AFM 
decimations. Since each decimation is equally likely to 
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be AFM and FM, the probability of having precisely k 
FM decimations is is 2~'''~^. The expected I is therefore 



<;>= 



OO 

E 

ra=0 



15 
12 



1 



15 
12 



(61) 



C. Entropy Calculation 

We now use the knowledge of mean decimation times to 
do the entropy calculation. We are interested in the scal- 
ing limit of large L, which translates to looking at large 
I = log Y- ■ In this case a complicated tree-like structure 
forms over each endpoint of the length L interval, and we 
need to figure out its entropy contribution. As we men- 
tioned above, for the case of the random Fibonacci chain 
with only AFM couplings, where the picture is a rain- 
bow diagram straddling each endpoint, the asymptotic 
contribution in the large rainbow limit of each singlet in 
the rainbow is log2 r. We will now find an analogous 
rainbow picture for the mixed Fibonacci chain. 

Let us focus on just one boundary of the interval, so 
we have one partition bond. The ground state trivalent 
graph that forms over it can be decomposed into con- 
nected tree components (fig. 11a) as follows. The first, 
innermost tree is generated by all the FM decimations 
prior to the first AFM decimation of the partition bond 
(if the first decimation is AFM, the tree is just a singlet). 
The next tree is generated by all the FM decimations 
between the first and second AFM decimations, and so 
on. As shown in (fig. 11a) these trees can be thought 
of as "thickened stripes" - it's just that now the stripes 
consist of not only a singlet, but an entire tree straddling 
the partition bond. The idea now is to use _F-matrix 
relations to decompose each tree into a superposition of 
graphs which have only or 1 r lines straddling the par- 
tition bond, so as to get the ground state to look like a 
superposition of rainbow diagrams (fig. lib). The en- 
tropy will then be a sum of a rainbow contribution and a 
contribution coming from the entanglement between the 
graphs on either side of the partition bond joined by the 
rainbow stripes. 

Before we go into details, let's compute a specific 
example. Consider the graph in fig. 12, which de- 
scribes two FM decimations followed by an AFM one. 
To compute its entanglement entropy we apply the F- 
matrix move as shown in fig. 12 and decompose it 
into a superposition of the two graphs on the right 
side of the equation. The entanglement entropy is then 



'log2 



'^log2(7 



In general we'll need to 



apply many i^-matrix moves and the superposition will 
be more complicated. 

Let's proceed step by step. Label the "thickened 
stripes" (trees) in fig. 11(a) by an index i, running 
from 1 to the number of thickened stripes n, and sup- 
pose the ith stripe connects a region Ai on the left side 
of the partition to a region Bi on the right side. Con- 
sider Hilbert spaces HAi and iJ^. associated with these 




Ai 



A3 A^ A, 



FIG. 11; a) Trivalent graph representing a ground state gen- 
erated by the real space RG. It can be decomposed into tree 
components, which are separated by AFM decimations of the 
partition bond, denoted in the figure by dotted lines. There 
could also be tree diagrams that don't straddle the parti- 
tion bond and hence do not contribute to the entanglement 
entropy; we omit them from the illustration for clarity, b) 
After applying F-matrix relations we can reduce the trivalent 
graph in a) to a superposition of graphs of this form. Here 
the dashed lines denote either a 1 (trivial) or r (nontrivial) 
line. 
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FIG. 12: Application of one _F-matrix move decomposes the 
graph on the left into the superposition on the right. 



sites. These are spanned by trivalent graphs having end- 
points on those sites, as before, but this time, because 
the regions may have nontrivial topological charge, we 
have the familiar decompositions Ha^ = FfAi'^ © -ff^i/ 
and Hsi — Hsi © -ffs; . The index and 1 just corre- 
sponds to whether or not Ai and Bi are connected by a 
T line. The Hilbert Space Hi of the union is 



Hi — Ha 



Hf 



Ha^<S>He 



(62) 



The ground state is a product state in H = ^- Hi . 
Using the decomposition for Hi we write the factors 
tpi = Oiiipi + Piipi ■ Here ipi'' are normalized to have 
norm 1 and \ai\'^ + |/3ip = 1. For convenience we take 
ai and (3i real and positive. Let 7 be the average over 
i of |/3ip. When we foil the above product we get the 
ground state as a superposition of 2" states with differ- 
ing rainbow configurations. One can check (by taking 
inner products and using the no tadpoles rule) that not 
only are all these states orthogonal, but all the states 
that enter into the Schmidt decomposition of one (on, 
say, the left side) are orthogonal to all the states that en- 
ter the Schmidt decomposition of the other. This yields 
a block diagonal decomposition of the density matrix of. 
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say, the left side. Thus we can deal with the blocks in 
this block diagonal decomposition separately. We label 
these 2" blocks with a label b. Each block b corresponds 
to a choice oi hi = 0, 1, where hi specifies the topological 
charge of Ai . The trace of such a block b is 



h^Y[, 



2hi Q2(l-hi) 



(63) 



Let's compute -Tr (Mf,logMf,) for block b. To do this 
we choose a convenient basis for the space where the, 
say, left components of the Schmidt decomposition of the 
state corresponding to this graph lie. In general, a ba- 
sis can be given by the set of all labelings of a trivalent 
tree, consistent with the fusion rules (i.e., you can't have 
two I's and a r at a vertex). So we pick a tree for each 
region Ai, and then join these up as in Fig. 13, for a 
tree defined over the whole left region (note that we're 
in the subspace where the topological charges of each Ai 
are fixed by hi, so we're not looking at all labelings, but 
fixing some of the edges to be r). This is precisely the 
kind of graph used in 21 to compute the entropy of the 
rainbow diagram. The only difference in our case is that 
we have extra degrees of freedom corresponding to the 
graphs for each Ai. So we've found a basis which con- 
sists of labelings of several subgraphs of a trivalent tree 
which do not interfere with each other. Namely, these 
subgraphs are the graphs near each Ai , and an extra one 
consisting of the edges which link the various Ai - it will 
lead to the "rainbow" contribution in the equation be- 
low. Because there are no mier-sub-graph constraints, 
i.e., the labelings on the sub-graphs can be chosen inde- 
pendently, the entropy is simply a sum of contributions 
from each subgraph. 



-Tr (MblogAfb 



-U log tb + h Sb 



raiiibo 



tb 2_^ ^i 



(64) 
Here the first term comes from the normalization of the 
block b, which has overall trace tb, the second term is the 
"rainbow" contribution mentioned above, asymptotically 
equal to log2 r times the number of nonzero hi, and Si-' is 
the entropy associated to ipi-' (i.e the contribution from 
the graphs around Ai and Bi. Summing over all blocks 
b and performing some elementary algebra, we get the 
entropy to be 

S = J2{-2a,^\oga, -2/3,^ log f3, + a,^S^°+f3,^S,^) 

i 

+ <56™"'^°"> (65) 

where the average < S^'^™'°°'" > is taken over all blocks 
b with weight tb. This average is approximately equal to 
log2 T times the average number of stripes, 771, so that 

S - ^(-2a,2loga, -2/3,2 log/?, +a,25,"+A'5,i) 

i 

-\- 7nlog2T (66) 
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FIG. 13: A trivalent graph whose labelings give a basis for 
the Hilbert space left of the partition bond. Here we have a 
tree for each region Ai, with the line leading out of Ai fixed 
to be r (for i = 1, 3, 4) or 1 (for i = 2). The rest of the lines, 
including the lines connecting different Ai, can be labeled at 
will, consistent with the fusion rules. The label of the dashed 
line gives the total topological charge of the left side of the 
system. 



Recall that we defined 7 above to be the average over 
i of |/3i| , which is just the fraction of regions Ai which 
have non-trivial topological charge. 

In eqn. 66 the first quantity sums up the tree contri- 
butions S'trec discussed at the beginning of this section. 
Thus eqn. 51 shows that it is equal to 



- < S'trec > Ini. 



(67) 



To compute < S'trec > we note that the average is taken 
over all possible trees generated by FM decimations be- 
tween two consecutive AFM decimations. There are 
many possible trees, since we can choose the number 
r of FM decimations, and for each FM decimation we 
must decide whether to decimate with the right or the 
left site. The probability of each such tree is 2"^'""^. We 
have computed this average numerically via a Mathemat- 
ica program. The program basically takes each possible 
tree and builds up the corresponding wavefunction "0 step 
by step in a convenient basis by applying the FM dec- 
imations. It then traces out half the system, finds the 
eigenvalues of the density matrix, computes the entan- 
glement entropy, and finally averages over the trees. We 
obtain < S'trec >= 0.115 approximately. The program 
also computes 7 = 0.927. Putting this into (66) we get 
S'= I InL (0.115 -H 0.927 log2(T)) so that 



S ■ 



0.234 log2 L. 



(68) 



The program goes up to r = 9, and we can bound 
the error obtained by omitting the remaining trees by a 
quantity exponentially small in r. Basically this is be- 
cause the probability of having a tree with a given value 
of r is exponentially small, 2''+^, whereas the maximal 
entropy contribution of such a tree only scales linearly in 
n (because the dimension of the Hilbert space is expo- 
nential in n). This argument yields a rigorous bound of 
±0.0006 on the coefficient in (68). 

Thus the effective central charge we obtain for the 
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abelian system: 




FIG. 14: a) A rainbow diagram, b) Typical state at the FM- 
AFM point. The dashed lines denote the location of AFM 
decimations, which separate the "thickened" stripes. 



mixed fixed point is: 



^eff 



3 • 0.234 = 0.702. 



(69) 



As we will discuss in the conclusion, this result is big- 
ger than the effective charge in the antiferromagnetic 
random-singlet fixed point. 



V. NON-ABELIAN ANYONS AND INFINITE 
DISORDER 

A. Fusion rules and real-space RG 

In this paper so far we analyzed the infinite disorder 
fixed points of the Fibonacci model. Rather surprisingly, 
the Fibonacci anyons lend themselves very readily to real- 
space RG analysis, and gives rise to a remarkably rich 
phase diagram. It is then natural to ask: will similar 
effects arise in other chains of non-abelian anyons? 

In fact, a deep relationship exists between real space 
RG and the behavior of non-abelian anyons. To see 
this consider first the ground state of conventional spin 
chains. In order to find the ground state of a conven- 
tional spin chain in a strong-disorder phase, we would 
apply real-space decimation rules to bond with strong 
coupling. The type of decimation we apply will depend 
on the local Hamiltonian and the Hilbert space of the 
system. For instance, in the spin-1/2 (AFM) Heisenberg 
model, two neighboring spin-l/2's can fuse according to 
the the SU(2) rule: 



i«i = 0®l. 



(70) 



A decimation rule applied to these two neighboring spins 
will choose one of the fusion subspaces - the spin-singlet 
or spin-triplet - according to the local Hamiltonian. The 
generalization of this principal to the case of non-abelian 
chains is nearly trivial. The spin-compounding rule, Eq. 
(70), is substituted by the fusion algebra of the non- 



a®b — 



'E 



n: 



ab'- 



(71) 



where iV^^ is the number of ways the superselection sec- 
tors a and b can fuse into c. 

A major difference, however, between rules (70) and 
(71) is that fusion rules for a non-abelian algebra are al- 
ways closed, while in regular spin-chains, the fusion rules 
include an infinite set of subspaces. The closure of the fu- 
sion rules for non-abelian anyons is a manifestation of the 
nonlocality of their Hilbert space, and therefore unique 
to these systems. It implies that one can always apply a 
real-space RG scheme without ever generating new types 
of coupling in the renormalized Hamiltonian. Further- 
more, just as in conventional spin chains, a decimation 
will result in either in a Ma-Dasgupta renormalization of 
the neighboring couplings, or in their multiplication by a 
factor of magnitude smaller than 1. Therefore sufficiently 
disordered (and most likely even weakly disordered) non- 
abelian chains will exhibit an infinite randomness behav- 
ior in the large length scale properties of their ground 
state. 



B. S > 1/2 Heisenberg chains and the SU{2)k fusion 
algebra 

The above observation is easily demonstrated using the 
mixed FM/AFM fixed point of the Fibonacci anyons. 
Both FM and AFM couplings between two Fibonacci 
anyons lead to fusion into either a Fibonacci anyon, or 
the vaccum: 

r (g) r = 1 ® r. (72) 

Therefore we can generically carry out a real-space RG 
analysis to its conclusion. But in spin-1/2 chains with 
nearest-neighbor couplings that could be either FM or 
AFM, it is easy to see that we generate higher and higher 
spins, and as a result do not flow to an infinite random- 
ness fixed point (although a fixed point of the mixed 
spin-1/2 chain was observed numerically in 24,25), un- 
less their Hamiltonian is restricted, e.g., by a symmetry 
in the problem which prevents large-moments formation. 
This is the case in 5 > 1/2 Hiesenberg models^^"^'* , which 
we will now briefly discuss. 

Disordered Heisenberg spin chains with spin S > 1/2 
were successfully analyzed by a real space decimation 
procedure that instead of forcing two strongly interact- 
ing sites into their lowest energy subspace (usually the 
singlet), just forbids them from their highest energy sub- 
space, (usually with spin 25'). This gives rise to sites 
becoming effectively lower-spin sites, with spins Si — 
1/2, . . . , 5^®'^^. Although the bare Hamiltonian contains 
only antiferromagnetic couplings, the decimation proce- 
dure also generates ferromagnetic bonds. These raise the 
spectre of large-spin moment formation, but the bipar- 
titeness of the chain in the bare Hamiltonian guarantees 
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that these FM couphng can never give rise to a spin larger 
than the original spin. 

The fact that in disordered Heisenberg models the real- 
space decimation rules only allow the formation of spins 
no larger than the original spins, makes these rules almost 
identical to the fusion rules of the truncated SU(2) rep- 
resentations, SU{2)k with k = 2S. As an example con- 
sider the spin-1 Heisenberg model. After some renormal- 
ization, the spin-1 chain effectively contains sites com- 
pletely decimated, spin-0, partially decimated, spin-1/2, 
and sites that are spin-1. Upon reals-space decimation, 
the ground state is formed by the following 'fusion' rules: 



101 = 

1(g) 1/2= 1/2 

1/201/2 = 00 1 



(73) 



These fusion channels are picked energetically; i.e., two 
spin-1 sites can fuse into a spin-2 moment, but this will 
be very costly, and will be excluded from the ground state 
wave function. The two fusion possibilities of the spin- 
1/2 in the last line indicate that spin-l/2's can have FM 
and AFM interactions. If we now compare this to the 

su2{2y. 



e(g)e=l 

e(g) cr = (T 

CT(g)Cr = l©e 



(74) 



we can identify the non-trivial superselection sector, ct, 
with spin-1/2, as expected form the Bratteli diagram, 
and the trivial sector e with spin-1. 

Indeed the two fusion rules are essentially identical. 
But as opposed to rules (73), which are imposed by en- 
ergy consideration, the fusion rules (74) are complete, 
and describe the full Hilbert state, rather than the 
ground state. Therefore the ground state of a disordered 
SU2{2) chain is different than that of a spin-1 disordered 
Heisenberg chain. In fact, the SU2(2) reduces to a ran- 
dom Majorana chain, analyzed in Ref. 21, which, has a 
random singlet ground state. A similar situation prevails 
in the case of the spin-3/2 Heisenberg model^*: the dec- 
imation rules for the spin model are almost exactly the 
same as the fusion rules for SUs{2) (except for the Heisen- 
berg model not allowing the fusion 101 = 1, which could 
be corrected by allowing biquadratic coupling). Never- 
theless, the spin-3/2 Heisenberg chain exhibits two ran- 
dom singlet phases, separated by a 4-domain permuta- 
tion symmetric fixed point, while the behavior of the 
non-trivial sector of the 31/3(2) is given by the above 
analysis of the Fibonacci chain. 



C. Novel infinite randomness universality classes in 
non-abelian anyons? 

Although the analogy between the spin 5* > 1/2 
Heisenberg model decimation rules, and the fusion rules 
of SU(2)k theories does not help us find new ground 



states, it demonstrates something rather important. Just 
as the random Heisenberg models allowed the discovery 
of the permutation symmetric sequence of infinite ran- 
domness fixed points^"'^, we expect that an investigation 
of disordered SU{2)k chains will also lead to novel infi- 
nite randomness universality classes. We leave this study 
for future research. 



D. Limit oi k ^ 00 

In the case of the disordered AFM golden chain, it is 
possible to generalize the setup slightly by considering 
other larger values of k. Retaining the quantum spin 1 
representation for the anyons out of which we build the 
chain, the entire analysis of 21 goes through, with only 
the quantum dimension changed from t to 2cos(7r/(A: -|- 
2)). It is gratifying to see that in the "classical" limit k -^ 
00, we reproduce the spin 1/2 result, with each singlet 
contributing logg 2 = 1 to the entanglement entropy. 



VI. CONCLUSIONS 

In this paper, we carry out an exhaustive analysis of 
the simplest random chain of non-abelian quasiparticles: 
the Fibonacci, or golden, chain. Using real-space RG, we 
are able to analyze the phase diagram and stability of the 
entire parameter range of the nearest-neighbor Fibonacci 
chain, where each pair of neighboring sites interacts by 
assigning an energy cost for fusing in the trivial channel 
or in the anyonic channel. 

The phase diagram we find is split between two phases, 
both of which are infinite randomness phases. When 
there are only couplings favoring fusion into the trivial 
channel (i.e., only AFM couplings) the flow is to the ran- 
dom singlet fixed point. When any finite density of 'ferro- 
magnetic couplings', i.e., couplings preferring the r fusion 
channel, are sprinkled in, the random singlet fixed point 
is destabilized, and the chain flows to a mixed infinite- 
randomness phase, which is characterized by the energy- 
length scaling exponent tp — ^, and a coupling distribu- 



tion function p{J) oc 



with X = 2, for both FM 



\j\i-x/ri 

and AFM cite. This stable fixed point, somewhat surpris- 
ingly, is in the same universality class as the transition 
point between the Haldane phase and the random singlet 
phase of the random spin-1 Heisenberg model^^'^^'^^. For 
the golden chain, the mixed fixed point also has a dia- 
grammatic representation in terms of random trivalent 
graphs. The mixed fixed point is the first non-singlet 
stable infinite randomness fixed point to be discovered. 

Another important character of this new IR fixed point 
is its entanglement entropy. Both infinite randomness 
fixed points exhibit the characteristic logi scaling: 



S 



1 



mixed 



InL. 



(75) 
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The coefScient in front of the In for the pure AFM chain 
was computed in Ref. 21, where it was found that it re- 
flects the quantum dimension of the Fibonacci anyons: 
Inr = 0.481. The entropy scahng calculation in 



respect to the FM one, unless another critcal point ap- 
pears in between, which we speculate is unlikely. On 



-eff 



the mixed phase is more intricate owing to the compli- 
cated trivalent graph nature of the ground state. We 
found the effective central charge, c™/?^'', to be: 



^myed ^ 3 . Q_234 = 0.702. 



(76) 



Since this result was obtained through a combination of 
numerical and analytical methods, it is hard to gain an 
intuitive understanding of the numerical result. Never- 
theless, it is interesting to compare it to its pure-system 
analog, and to the effective central charge of the ran- 
dom singlet phase. It is most likely that the mixed IR 
phase is also the terminus of flow from the ferromag- 
netic pure fibonacci chain. The central charge of the 
critical FM golden chains was determined in Ref. 9 to 
be c = 4/5 = 0.8 > djff"^. Hence the effective cen- 
tral charged dropped along the flow. Comparing our re- 
sult, though, to the central charge in the random sin- 
glet phase immediately reveals that the effective central 
charge increased in the strong-randomness RG flow from 
the random singlet phase to the mixed IR phase. Thus 
the suggestion that strong-randomness flows may have a 
c-theorem associated with them is contradicted. 

This result is rather novel, and it is worth mentioning 
that one can see that it is true without having to do the 
full calculation of the entanglement entropy in the mixed 
FM/AFM phase. To see this, one first of all notices that 
the average < I > between AFM decimations, equal to 
3/2 (eqn. 61), is half of that in the AFM phase. This 
means the treelike structures form twice as fast as the 
singlets in the AFM phase. Now, to compare coefficients 
in front of In L one must multiply by the energy length 
scaling exponent, which is 1/3 in the mixed phase and 
1/2 in the AFM phase. So the number of treelike struc- 
tures forming over the endpoints is 4/3 times the number 
of singlets in the forming in the AFM phase. However, 
with probability 1/2 the treelike structure is simply a 
singlet (AFM decimation following another AFM deci- 
mation), and with probability 1/4 it is simply a tree on 
3 sites (one FM decimation between the AFM decima- 
tions), whose contribution to the entropy is the same 
as a singlets. Thus the treelike structures contribute at 
least 3/4 as much entanglement entropy as singlets, and 
given that there are 4/3 times as many of them as sin- 
glets in the AFM phase, the mixed phase has at least 
as much entropy as the AFM phase. The fact that the 
more complicated trees have a nonzero contribution im- 
mediately shows that the entropy of the mixed phase is 
in fact higher than that of the AFM phase. 

Now, when we contrast the central charges of the of 
the two critical phases of the pure chain, c^^^ = 0.7, 
and c^^ = 0.8, we find that by the Zamolodchikov c- 
theorem,^^ the AFM phase must be a stable phase with 



Disorder 



:0.702 




FM/AFM 



FIG. 15: Flow diagram of the pure and disordered golden 
chain. In the pure chain, assuming no intervening fixed points 
exist, the FM fixed point is unstable to flow to the AFM fixed 
point, as inferred from the Zamolodchikov c-theorem. In the 
disordered chain, however, the fiow is in the opposite direc- 
tion, with the mixed FM/AFM phase, which is most likely 
the terminus of the fiow from the pure FM phase, being sta- 
ble relative to the random singlet phase, which is the result of 
disordering the pure AFM phase. The fixed point (effective) 
central charges are also quoted. 



the other hand, the flow in the random golden chain is 
the opposite: the mixed FM/AFM phase is stable for es- 
sentially all chain coupling distributions, except for the 
point in which all couplings are antiferromagnetic. This 
situation is summarized in Fig. 15. 

Most importantly, we also observed in this paper the 
close connection between a fusion algebra and real-space 
RG. This connection implies that essentially all strongly 
disordered phases of non-abelian chains will be of the 
infinite randomness class. In the future we intend to 
analyze the random phases of non-abelian chains with 
different fusion algebras. While this research is intended 
to expand our understanding of random non-abelian sys- 
tems, it may also lead to the discovery of new infinite- 
randomness phases and universality classes, beyond the 
permutation-symmetric sequence of Damle and Huse^^. 
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